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Abstract: We report on the computation of a class of massless bosonic three-loop vacuum 
sum-integrals which are key building blocks for an evaluation of the Debye screening mass in 
hot QCD. Generalizing known techniques and introducing the concept of tensor reduction by 
dimensionality shifts (known to the zero-temperature community since the work of Tarasov 
in 1996) to finite temperature, we are able to treat hitherto unaccessible cases, which will 
allow us to finalize the long-term project of NNLO Debye mass evaluation. 
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1 Introduction 

Taking as a motivation the startling imbalance between highly developed analytic, systematic 
algorithmic as well as numeric methods for multi-loop integrals in zero-temperature field 
theory on the one hand (see, e.g. [1]), and the many unsolved technical challenges one is 
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confronted with when dealing with phenomenological problems in cosmological or heavy-ion- 
collision related contexts that are formulated within finite temperature field theory on the 
other hand (see, e.g. [2-6]), we continue our line of systematic investigation of the latter class of 
problems [7-12]. Once more diving into a complicated and rather technical issue, our aim is to 
exhibit and develop generic tools that allow for progress on the thermodynamic (equilibrium) 
front. Nevertheless, while in this note we merely focus on a specific class of multi-loop 
sum-integral that arise in the thermodynamic setting, our final result will find a concrete 
application in the determination of matching parameters (such as the Debye screening mass) in 
the systematic program of an effective theory treatment of hot QCD thermodynamics [13, 14]. 
We relegate the finalization of the long-term project of determining the Debye screening mass 
of hot Yang-Mills theory to NNLO [9] to an upcoming publication [15] , and here focus 
on the well-separated problem of evaluating the last missing sum-integral in dimensional 
regularization, up to its constant term. Let us just mention here that, once full 3-loop results 
for m\ (interpreted as a matching coefficient between 4d thermal QCD and 3d electrostatic 
QCD (EQCD) in the framework of dimensionally reduced effective theories [13, 14]) are 
available, a certain contribution of order g to the hot QCD pressure can be deduced, and 
refer to Sec. 6 of Ref. [9] for a more detailed discussion of this potential phenomenological 
application. 

There are a number of well-established techniques that prove useful when evaluating 
sum-integrals at higher loop order. Starting from the seminal work of Arnold and Zhai [2], 
a number of specific cases have appeared in several applications, ranging from 3-loop sum- 
integrals needed for QCD thermodynamics [2, 3], to even some 4- loop cases that have served to 
quantify high-order corrections in the thermodynamics of scalar theory [16, 17] as well as large 
Nf QCD [18]. In each case, the corresponding sum-integrals that appear in the calculations 
have been evaluated to sufficient depth in their e-expansions, building a (small but valuable) 
database of master sum- integrals. There is an urgent need to enlarge this database, e.g. in 
order to make progress on the precision of matching computations [9] (involving moments of 
3-loop self-energies) or on the determination of the physical leading order (for a justification 
of this term see e.g. Sec. 6 of [9]) of hot QCD observables [4, 10] (where further 4-loop masters 
are required). 

More recently, in a series of works we have re-examined these computational techniques 
as well as most of the known specific cases of sum-integrals, and put them on a more generic 
footing, using notation that allows for generalizations [8, 9, 11, 12]. Performing those gener- 
alizations for computing previously unknown cases [12], we have pointed out that integration- 
by-parts (IBP) relations can be used profitably here as well. IBP relations have already been 
successfully applied in the thermal context to the reduction process itself [9, 10], but now 
they can also be employed to help in evaluating the corresponding master sum-integrals that 
remain after the reduction algorithm has halted - mainly in order to transform infrared di- 
vergent parts of dimensionally regularized masters in terms of convergent ones [12]. We will 
see that also in the present context, they play an important role. 

In this note, we wish to infuse another new idea to the art of evaluating sum- integrals, 
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namely to use the ideas of Tarasov [19] for tensor reduction. The main goal is to avoid the 
traditional approach of projection methods that - if applied to thermal field theory where the 
rest frame of the heat bath breaks rotational invariance and introduces a vector U = (1, (0)) 
that all tensors can depend on - lead to inverse structures (such as 1/p 2 ) of a different form 
than propagators, hence leading outside the class of sum-integrals that one has started with 
(in fact even outside its natural generalization as suggested by IBP [10, 11]). The price to 
pay is a shift (in units of two) of the dimensionality of the integral measure, and an increase 
of powers of the propagators that are present in the corresponding sum-integral [19]. We 
will argue that this price is not too high, and show explicitly that one can deal with these 
dimensionally shifted cases. 

As a concrete example on which to test our new methods, we will consider a specific mass- 
less bosonic three-loop sum-integral of mass dimension two (which we shall name ^3,-2) here, 
which constitutes one of the last building blocks needed for evaluating the Debye screening 
mass of thermal Yang-Mills theory. In fact, in order to stress the generality of our method, 
we will treat a more general (infinite) class of sum-integrals (Mn,-2, f° r integer values of N) 
for most part of the paper, of which ^3,-2 is just a special case. 

The paper is organized as follows. In Sec. 2, we introduce the class of 3-loop tensor sum 
integrals Mn,-2 on which we wish to test our new idea of tensor reduction, and explain its 
decomposition into scalar pieces. Owing to the structure of the latter (they all contain two 
one-loop sub-integrals and are hence of so-called spectacles-type) we further decompose them 
into pieces that either allow for analytic solutions, or are finite such that they can be evaluated 
numerically. Sections 3 and 4 then deal with the evaluation of non-zero (Matsubara) modes 
and zero-modes in turn, specializing to the specific sum-integral M3-2 for concrete results, 
which are then summarized in Sec. 5. We conclude in Sec. 6 and have relegated some technical 
material to the appendices. 

2 Decomposition of M. n-2 

The sum-integrals M j\r -2 that we shall be concerned with here can be represented as a subset 
of a more general class as defined by (see Eqs. (A. 23), (A. 30) of Ref. [3] and the review [11]) 



We employ (Euclidean, bosonic) four-momenta P = (Pq, p) = (27rra p T, p) with P 2 = Pq +p 2 ; 
the temperature of the thermal system is T; and the sum-integral symbol stands for 



The known non-trivial instances of the class Eq. (2.1) are A^o,o an d M2-2 [2] of (mass-) 
dimension four, as well as the dimension two case M\q [20], all of which have been re- 
evaluated in unified notation in [11]. Furthermore, as has already been pointed out in the 






with d = 3 — 2e . 



(2.2) 
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latter reference, the numerators of .Mat-i can be eliminated by using the denominator's 
invariance under momentum shifts Q — > —P — Q and R — > — P — R, resulting in 

Mn --' = t Q IPT v\p + Q? ~\ M "- l ° ' <2 - 3) 

where the first term on the right-hand side (rhs) has factorized into a trivial 1-loop tadpole I\ 
(cf. Eq. (A. 2)) times a two-loop sum-integral of sunset-type which can be further factorized 
into a product of one-loop cases by IBP (see e.g. App. B) and is hence trivial as well, while 
the second term on the rhs is again in the class of Eq. (2.1), although without scalar products 
in the numerator and hence much simpler. 

The massless bosonic 3-loop vacuum sum-integral A4 j\r -2 is therefore defined as 

Mn _ 2S 4 h [M-*)T . (2 . 4) 

As mentioned above, one representative of this class, A^2,-2, is already known [2, 3], since 
it has entered the 3-loop result for the thermodynamic pressure of hot QCD. However, its 
determination (as reviewed in [11]) was somewhat contrived, mainly due to difficulties in 
treating its numerator structure - a problem that has triggered much of what we report here, 
and that we will alleviate Sec. 2.2 below. 

2.1 Splitting off trivial parts 

Expanding the numerator of the class of massless bosonic 3-loop vacuum sum-integrals AAn,-2 
inEq. (2.4) as [{Q-R) 2 ] 2 = A(QR) 2 -A(QR)(Q 2 + R 2 ) + Q A + R 4 + 2Q 2 R 2 , its spectacles-type 
structure becomes explicit when using the {scalar, vector, tensor}-like 2-point functions 

{iMP),n W P),iC(P) } s f Q jgjggg*, ■ (2.5) 

such that (omitting the argument (P) of all functions II for brevity from now on) 

Mn,-2 = £ jpip Wrn^ - 'ffl&nft - 4n^n^ + n_ u nn + n n n_ n + 2n 01 n i} 

7 -J w {4n^n^ - 2P 2 n n n 10 + 2n 10 n 10 } (2.6) 
p \ 



= ^ p |pip n n n n " 2J^ j^yhl + 2 Wi u , (2-7) 

where in the second and third line we have used some trivial properties of the scalar and 
vector-like 2-point-functions, which derive from the shifts Q — > —P — Q or Q — > —Q and read 

ru = u ba , n a0 = i° a , noo = ig = i l = o , n_ la = p 2 i° a + i° a ^ ; (2.8) 

Tq 

Kb = -P'Uta - Ka . = ~ ^aa , K = , Ka = 4° • (2-9) 
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In Eq. (2.7), the third term is a product of 1-loop sum-integrals and hence known analytically, 
cf. Eq. (A. 2), the second term is a 2-loop problem and hence equally trivial (it factorizes into 
1-loop integrals via IBP, see Eq. (B.6)), while the first term needs further treatment and will 
be the focus of the next section. 



2.2 Taming the tensors 

To proceed, rewrite n^Tl^ 



n oo n oo , on * TT 0i 4- J] ij Tl ij 
J- 1 !! 11 !! + Zli ll li ll + 1L n LL n 



L1 . Noting that standard tensor 
reduction of the 3-vectors present in the numerator of e.g. = |lj + ^ IIn — 2Il5] ) |- 
produces inverses such as 1/p 2 , which we want to avoid since they are not contained in 
the structure of the original sum-integrals we started with, we choose to employ Tarasov's 
T-operator technique [19] to trade scalar products of 3-vectors for higher dimensions. 

To this end, we regard the (massless, 4d) sum-integrals as sums over corresponding 
(massive, 3d) spatial integrals, taken at specific values of the masses. In the spirit of Ref. [19], 
it is then advantageous to generate the irreducible scalar products of 3-vectors (qr in our case) 
by derivatives acting on a generating function exp(— 2aqr): 



1 



i 

[P2p 



T-rOi ttOi 



if- 7 ' u ij 



T 3 ^2 Qo Ro d~2a ^nii (-Po, Qo, Ro, Po + Qo, Po + Ro', 

PaQoRo 



a=0 



T 3 9- 
PoQoRo 



2a 



^ivnii (-^0' Qoi R°i p° + ^5°' p° + R°> a > 



a=0 



(2.10) 
(2.11) 



with the generic (massive, 3d) 3-loop vacuum integral 

-(d) e -2a(qr) 



pqr 



-ml] u i [q 2 +m|] !y 2 [r 2 + m|] v 3 [( p +q)2 +m 2]i, 4 [(p + r .)2 +m 2]„ 5 



This integral has another useful representation, whose a-derivatives will become simple: intro 
ducing Feynman parameters for the propagators 1/A U = f °°d 



a u 1 e aA /T(v), completing 



squares in the exponential, shifting momenta and rescaling them, it follows that 

/ 5 f°° a Ui ~ x e _Qim ? \ r 

' da,; 1 — ; I D(aj, 



Ad) 

± Vl...5 



n 

vi=l 







D(oij) + 020405 — a (01+04+05" 



where D(a,j) is the graph polynomial [21] 

D(ctj) = (o 2 + 04)0305 + (03 + 05)0204 + oi(o 2 + 04X03 + 05) 



-d/2 



(2.12) 



(2.13) 



The Gauss-integral in Eq. (2.12) depends on the chosen integral measure; in ours, it reads 
I 



d d p 
(2w) d 



= (47r) 3d / 2 . Now, note that a-derivatives acting on Eq. (2.12) lower the 
power of the last term by an integer n (which can be interpreted as a shift d — > d + 2n), 
while producing additional polynomial pre- factors in the Feynman parameters Oj . The latter 



- 5 - 



can then be pulled out of the integral letting ay — > d_ m 2 , and finally be absorbed in positive 
shifts of the propagator powers fj. We need the two instances (d = 3 — 2e) 



u - 2al Nllll 



0: 



2a J 7Vllll 



o=0 



o=0 



d 



,2(9 



u -'Afllll 



o=0 



d 75+ rW 

2 J iV1122 



o=0 



<9_ m 2+d_„,2+<9_ 



Willi 



o=0 



2 9! 



-n++r( rf ) 



7V1111 



o=0 



^ J Afll21 



^ J Aflll2 



o=0 



+ d(d + 2) £> ++ J 



7V1133 



o=0 



where V+I^ = (An) 3 1^ 2 ^ etc., such that Eqs. (2.10), (2.11) finally read 

[p2]7V 



P 



[P 2 } N 
1 



T-rOi rrOi 

ii ll li ll 



il ll li ll 



2 



— — — — + - T> + 

[p2]7V+l ^ 2 



jp2l7V —ii—ii 

and we obtain our final representation of Eq. (2.7) as 



^ + d(d + 2)V++ 



[p 2 Y 



n 3 in 31 

[p2]7V 



M 



N,-2 



n??n?? 



p p2N +4^^ + A^+ 



n n n 



ii 



p27V+2 



+ 2d£>~ 



n 21 n 



ii 



p2N 



+ 4d(d + 2)£> H 



n 3 in 3 i 



-3U*31 „,Ovf JJu ,orOrO r O 
p2lv 27 14- p2JV-2 + ll Nhh 



rO r O r O 



+ 



(2.14) 



= W 3 + MV A + NdV 5 + 2dV 6 + Ad{d + 2)V 7 - 21^ + 21% . 

It remains to calculate 1/3/1 3- loop spectacles- type sum-integrals {V^, V4, V5, Vq, Vj} in d/d+ 
2/d+A dimensions. Note however that their structure is quite similar, such that we will be 
able to employ a quite generic strategy for their evaluation, as will be explained below. 

2.3 Decomposition of spectacles-type sum-integrals 

Any spectacles- type sum- integral with two generic 2-point functions Ili(P), n 2 (-P) can be 
identically rewritten as (suppressing indices of IT; for a full definition, see Eq. (3.13)) 

fp [P 2 Y 



v = 
v l = 
v d = 
v z = 



nxn 2 
1 



P [P 2 ] n 

, (Po) m 

p [P 2 ] n 

Spo (Po 

p 



yf(inite) _|_ ^d(ivergent) _j_ -^z(ero) 

^ (Hi - nf )(n 2 - nf ) + (n x - nf )(nf - nf )} + {1^2}, 
n x nf - nf nf + (nf - nf )(nf - nf ) + \ nf nf } + {102} 



n,n 



[p2 ]„ "1 1X 2 [p2] 

where the n B (T = part; leading UV div), n D (integer power of P 2 such that [U B 
finite as e — )■ 0) and H c (sub-leading UV divergence and sometimes further subtractions) are 
chosen such that the terms in V { are finite (and could be evaluated numerically in coordinate 
space at e = 0); those in V d are simple 2-loop (first term) and trivial 1-loop structures; while 
V z is the zero mode treated in Sec. 4. The subtraction terms Y[ B ' C ^ D are detailed in Sec. 3.1. 



■ILII 



1^2 



(2.15) 
(2.16) 
(2.17) 
(2.18) 



n 



IS 
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3 Non-zero modes 



The goal of this section is to compute the V d of Eq. (2.17) analytically, and to evaluate the 
V f of Eq. (2.16) numerically, from simple low-dimensional integral representations. 



3.1 Generic formulae for 2-point subtraction terms 

For generic indices, let us specify the subtraction terms as 
1 g 00 (a,b,d + l) 



n 



'q [Q 2 ] a [{P + Q?\ h (p2)a+Md+l)/2 



with goo{a, b, d) = G(a, b, d) 



nf 61 = Up 



U^{P^A(P 2 )}=P Q A(P 2 ) 



Q [q^[(p + gn 

_ Po9io(a,b,d + l) 

(p2^a+b-(d+l)/2 

with 5io(a, b, d) = — G(a, b — 1, d) — G(a — 1, b, d) — G(a, b, d) 



(3.1) 

(3.2) 
(3.3) 



n 



ab2 = 



U,,U V 



q [Q 2 ] a [(P + Q) 2 ] b 



U^U U {g^A(P 2 ) + P^P U B(P 2 )} = A(P 2 ) + P 2 B(P 2 



g 2 i(a,b,d + l) Pq g 2 o(a,b,d + 1) 



+ 



(p2^a+b~(d+3)/2 (p2^a+b-(d+l)/2 

with 021 (a, M) = \2G(a-l,b,d) + 2G(a,b- l,d) + 2G(a- 1,6- l,d)- 

1 



G(a, b, d) - G(a - 2, b, d) - G(a, b-2,d) 



4(d-l) 



and g 2 o(a, b, d) = G(a — 1, b, d) — dg2i(a, b, d) 



+ 



1 b 



T\C — TTt 

(p2)b ' (p2^a ' abl — xl abl (p2y ' 



T 2 t2 i p2 T 

ttC _ ttB i fo , 1 b "+~ M) 1 b 
Ll ab2 = Li ab2 TT^TT + 



n 



D 

abc 



2\e 



(aT 2 )' 



(P 2 ) b (P 2 ) a ' 
nf fec (note that this is oc (^integer ^ for any ^ = odd - 2e) 



(3.4) 

(3.5) 
(3.6) 
(3.7) 

(3.8) 

(3.9) 



where the coefficient functions gij derive from 4d rotational invariance. Since efficient com- 
putation needs compact notation, we summarize the various II s and Tl c as given above as 



[c/2] 



[P ) c - 2n gc,n(^b,d+l) 



n=0 



(p2^a+b~(d+l)/ 



2-n 



T\C _ TT 

abc abc 



■ I + l " ir ' :( 1 + (-.rE( 2 ;)fw (pr! " 

n=0 ^ ' 



2 (P 2 ) 

where [c/2] = {c/2, (c — l)/2} for c = {even, odd}. 



{P 2 ) a ' 



(3.10) 
(3.11) 
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3.2 Generic formulae for divergences V d of non-zero modes 

Due to the structure of n^ c , cf. Eq. (3.9), Eq. (2.17) factorizes as 

+ {1^2}, (3.12) 

where the a.% (cf. Eq. (3.9)) are constants that might be chosen such as to simplify the finite 
parts V i , and which also serve to facilitate comparison with other approaches (e.g. [11, 12, 20]). 
For the generic spectacles-type 3-loop sum-integral 



^(si... 5 ;s 6 ...8) = 



[Po) S6 (QoY 7 (R 



I* 8 



tpQR [P 2 ] S1 IQ 2 Y 2 im s * l(P + Q) 2 Y* i(P + R) 2 ] ss 

= f p p£ n s2S4S7 (P) U S3S5SS (P) (3.13) 

we therefore write a general expression for the divergent parts of their non-zero modes: 

[«a/2] 

V d (si...8) = 9s s ,n( S 3i S5, d+1) f a2 (s 135 - ^ - n, s 2 , s 4 , s 68 - 2n, s 7 ; d) + 

n=0 

[»r/2] 

+ E 9s 7 ,n( s 2, »4, d+1) / ai (si24 - ^ - ra, s 3 , s 5 , s 67 - 2n, s 8 ; d) , (3.14) 

n=0 

(aj-independence follows when adding Vf T , cf. Eq. (D.35)) where the functions / are given by 

/«K. 5 ^)-^Pyr{-eEq. (3.12)} s23s 
1 



L'( Sl - e, S2 ...5, d) + /| 2 5 (d) (d, o)+ 



(aT 2 ) 

[«6/2] 

+ 

n=0 

1 

2 ys 5 ,n^2, 03, u, -r j.; ^ Sl23 _( d+1 )/ 



[«/2] , v 



+ ^fe,n(5 2 ,S3,d+ l)Cti 5 (d+l)/2-n( d »«) ( 3 " 15 ) 



with the one- and two-loop structures 

I b a (d, a) = I h a {d) - (aT 2 yn h a _ e {d) , (3.16) 
P a (d, a) = I b a (d) - 2(aT 2 )-% b _ e (d) , (3.17) 
L'( Sl ...5, d) = L(*i... 5 , d) - 5 S4 ^(s L ..3, a 5 , d) • (3.18) 

The latter reduces to 1-loop sum-integrals as L —>/■/, as can be derived systematically via 
IBP; the cases relevant for the present computation are given in appendix B. 
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3.3 Specific results for divergences V d of non-zero modes 

From Eq. (3.14), we obtain the desired expansions for the divergent pieces of {V3, V4, V5, Vq, V7} 
needed for the decomposition Eq. (2.14) of the specific sum-integral around d = 3 — 2e, 

T 2 (47rT 2 )- 3e 1 
e 2 288 



V£ = y a (31111;022) 
V£ = £> + V d (31122;011) 

V d = p+y d (4ini) 
v£ = p+y d (32in) 

V 7 d = V ++ V d (33311) 



(4tt) 



(4tt) 



73 



1+ ( Y§+7E + 241nG) e + 0{e 2 ' 



(4ttT 2 ) 



(4itT 2 ) 



-3e r 



-3e 



1 



(4vr) 4 

T 2 (AttT 2 )-^ - 
' (4vr) 4 e 2 96 
T 2 (4^T 2 )" 3e 1 

7 2 432 L 



( — — II7E — 120 In G + 24 ln(2vr) 



(4vr 



1 + 
1 + 



139 
l¥ 
26 13 



17 40 
y 7E -561nG+yln(27r; 



96 



■7E 



InG 



36 



ln(27r) 



where y d (41111) = y d (41111, 000) etc., the Glaisher constant G appears inside a logarithm, 
with 121n(G) = l + l)/£(— 1), and we refrain from listing the somewhat lengthy constant 
pieces, which contain the dependence on the arbitrary constants a% as defined in Eq. (3.12) as 
well. As mentioned above, this dependence will precisely cancel that in V { , which is contained 
in V*j, cf. Eq. (D.35). 

3.4 Finite parts V { of non-zero modes 

With the structure of the subtracted terms as in Eq. (2.16), they contribute to M.3-2 as 

M™'1 2 = 4V 3 f (31111, 022) + l2V + Vl (31122, 011) + 9X>+F 5 f (41111)+ 
+ 6P + y 6 f (32111) + 60P ++ F 7 f (33311) + 0(e) 
T 2 



[ni + 0(e)] 



m ps +0.0645513(1) , 



(3.19) 



where we discuss the numerical evaluation of m, for which coordinate-space methods prove 
useful, in App. D.l. The value given above was obtained setting all on = 167r 2 e 3 / 2-7 . 



4 Zero modes 

Here, we treat the Po = modes of all 3-loop sum-integrals of Eq. (2.14) that are needed for 
the specific case .M3 _2, while discussing the generic strategy. In terms of 



n. 



n 



and using the IBP relations of App. C to improve their IR behavior, they read 



Q 







216 



(4.1) 



V£ = 5(31111; 22) 



5(12121; 22)+5(12211; 22) - if A(221; 2, d) 
d-8 L {d- 5) 



-7 2 2 A(311;2,d) 



(4.2) 
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VI = £> + S(31122; 11) = after summation in <5p II a6 i , 
V£ = X>+5(41111) 

1/ 6 Z = X>+S(32111) 



d- 8 
1 



V 6 z -V + (l$A(All;0,d)) 



£>+S(22121) + ^2^+5(12221) - V+ (if A(321; 0, d)) 



d-5 L 

(d 2 - 13d + 38)£>++S(23321) - (d - 10)D++S(23222) 



(4.3) 
(4.4) 

(4.5) 
(4.6) 



^ = ^(33311) = 2 (d _ 2)(d _ 5 )(d-8) 
A generic decomposition into finite part (first line) and remainder of zero-mode integrals S is 



S(si... 5 ;s 6 ,s 7 ) = [n 



p [P 2 r 



[n a 



+ 



+ 



Spo n n A + 



[P 2 ] s l ± ^ S 2S4S6 ±± S3SsS 7 



[P 



Sp ° nd.._n. 



21si S 3S5S7 



where the subtraction terms are defined as 



gf(inite) _|_ ^d(ivcrgent) 



(4.7) 



U A bc = 6{d+2-2a-2b+c) 5 Po U B bc + 9(2a+2b-c-2-d) 6 Po "£ 

Tq 

= 6{d+2-2a-2b+c) S Po U B bc + 9{2a+2b-c-2-d) 5 Po 5, 



5 Q Qo 



TG(a,b,d) 



(4.8) 



■ (p2^a+b-d/2 ' 

which represents subtraction of the usual leading UV divergence H B given by the T = 
piece, whose analytical representation was given in Eq. (3.10), as well as subtraction of the 
zero- mode inside two-point functions II^o- 

4.1 Generic results for divergences 5 d of zero modes 

Noting that the subtraction terms n^ c are proportional to powers of p 2 , the last term of 
Eq. (4.7) vanishes identically since it is scale- free, while the other two can be expressed 
analytically using the 2-loop function A: 

5 d (si...5; s 6 , sr, d) = a(si... 7 , d) + a(si 325 476, d) + sca M re e , (4.9) 
a(si...7, d) = 6(d+2-2s 2 4 + s 6 ) e S6 g S6 , S6 / 2 (s2, «4, d + l)A(s 12 4 - d+1 ^ S6 ,s 3 , s 5 , s 7 , d)+ 
+ 8(2s24-se-2-d)5 SG TG(s 2 ,S4,d)A{si24 - ^,s 3 ,s 5 ,s 7 ,d) , (4.10) 

where e s = (1 + (— l) s )/2 is 1(0) for even(odd) index. 

4.2 Specific results for divergences V z,d of zero modes 



Collecting and expanding around d = 3 — 2e dimensions, we finally get for the divergent pieces 
of the zero modes 

ld T 2 (4vrT 2 )- 3e -1 



(4vr) 4 e 2 



48 



+ le + + io " 37E + 61n(27r) ' 



(4.11) 
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Va 



4 -0, 
d 



V, 



z,d 



rp2 


(4ttT 2 ) 


-3e _l 


1 + ( 


'5 


(4tt) 4 


e 2 


72 


v3 


rp2 


(47rT 2 ) 


-3e 5 


1 + 


f- 


(4tt) 4 


e 2 


144 






(47rT 2 ) 




1 + 


f- 


(4vr) 4 


e 2 


240 


Vl5 



7E + 41n(27r)J e + 0(e 2 
- 7E + 41n(27r) ) e + 0(e 2 ) 



(4.12) 
(4.13) 

(4.14) 

(4.15) 



Here, we have of course assumed that all six subtracted sum-integrals shown in the first line 
of Eq. (4.7) are finite. We prove this assertion by explicit computation in App. D.2, with 
results given in Eq. (4.16) below. 

4.3 Finite parts V z ' f of zero modes 

We now have to write down integral representations for the Il^-subtracted parts of our zero- 
modes and solve them numerically, to prove that they are finite indeed. The first line of 
Eq. (4.7), to be evaluated numerically at d = 3, contributes to as 



M 



4 

2 ~ 5 
6 

~ 5 
+ V~ 



f p 5 -0 { [n 222 - n£ 2 ] [n 112 - nf 12 ] + [n 212 - nf 12 ] 2 }- 
v + ^ p { [n 220 - i4 ] [n 110 - uf 10 ) + 2P 2 [n 220 - r4 ] [n : 



+ 



^{96 [n 320 - 14 ] 



[n 



310 



n 



310. 



+ 84 [n 



320 



n 



320. 



[n 



22 



(4vr)^ 



[n 2 + 0(e)\ 



n 2 » +0.24983747686(1) 



n&>] } 
n&>] } 

(4.16) 



where we discuss the necessary integral representations as well as their numerical evaluation 
leading to n 2 in App. D.2. 



5 Result for 

According to Eqs. (2.15)-(2.18) and Eqs. (4.2)-(4.7), each scalar spectacle-type sum-integral 
is decomposed as V = V { + V d + F z ' f + V z,d , and according to Eq. (2.14) (with Eq. (B.6) 
for the 2-loop case therein), the master integral A^3 5 _ 2 can be written exclusively in terms of 
such scalar spectacle-types and trivial 1-loop integrals. Collecting and expanding, we finally 
obtain 



with m 



T 2 


( V £ 1 


(4vr) 4 


\4irT 2 J e 2 


59 


83 2 


405 + 


T^lE +1E 


757 




7T 


2 --ln 2 2- 


2160 


3 



5 

36 ~ 
559 



(5.1) 



+ - ln(2vr) - - In G + — — t(3) 
1080 3 v ; 3 1080 sv ; 



1 



3780 



C(5) 



65 7 2 28 

— In G + — ln(2vr) - - ln(4vr) In vr + — 71 + 
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37 113 7 1 

+ — C(3) + — C(5) - — C'(3) + C'(5) + 2C"(-1) + ni + n 2 (5.2) 

3240 sv ; 226800 sv ; 540 s v 7 1890 s v; ^ sv ; t i t 2 v > 

« m + n 2 - 6.5129800346005775127 » -6.1985913(1) , (5.3) 

where the Glaisher constant G appears inside a logarithm, with 121n(G) = l + l)/£(— 1), 
and various zeta values as well as the Stieltjes constant 71, defined by £(1 + e) = 1/e + 7e — 
71 e + 0(e 2 ), enter the constant term. 

While it might at first sight seem that we have used an unnecessarily generic notation in 
our decomposition and treatment of the various contributing terms, let us remark that this 
was done in order to provide a setup that allows for solving many more sum-integrals than 
just the special case M.%-2 given above. Indeed, as we quickly demonstrate in App. E, the 
two previously known sum-integrals Mi t o [17] and V2 [12] (which each had required quite 
some effort to evaluate) are reproduced by our generic formulae, and there is no doubt that 
the same could be said about M.2,-2 [2, 3, 11], although we have not explicitly checked that. 

6 Conclusions 

In this work, we have made important progress in two respects. First, by transferring proven 
technology from zero-temperature field theory to the finite-temperature case (the T-operators 
of Sec. 2.2), we were able to map tensor sum-integrals onto scalar ones. Naturally, although 
we have tested this method on a specific case only, it is applicable much more generally. Much 
in the spirit of Ref. [22], by using the new tensor technique we have dissected the problem 
to its clean core, finding a class of spectacles-type sum-integrals that are amenable to known 
techniques. 

Second, we have worked out generic solutions for these massless bosonic 3-loop spectacles- 
type sum-integrals, which are applicable to a wide class of such integrals. Let us note that 
this is only the third instance of treating a more generic class of sum-integrals in a single 
computation (the first two such instance were recorded in [8, 9]), and represents the first steps 
beyond the case-by-case analyses found in the literature - a development urgently needed in 
order to close the glaring gap between well-established generic integration techniques in zero- 
temperature field theory and the few painstakingly derived cases that are known at finite 
temperatures. 

A generalization of our results to sum-integrals involving fermions should be possible (at 
least for vanishing masses and zero chemical potentials) in a straightforward manner, and 
could in fact even turn out to be structurally simpler due to the absence of zero-modes in 
fermionic lines, which had cost us quite some effort in the bosonic treated in Sec. 4 

above. There are a number of well-defined applications that await the evaluation of such 
fermionic sum- integrals, such as the iVf 7^ parts of EQCD matching coefficients as listed in 
[9]. 

As an application, of our generic formulae, we have given results for the new massless 
3-loop sum-integral M.3-2, which represents a key building block of the bosonic part of the 
Debye screening mass in hot QCD [9, 15]. Reassuringly, our generic formalism reproduces 
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the two previously known sum- integrals of the spectacles- type class without effort, as we have 
demonstrated in App. E. 
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A Standard integrals 

For convenience, we collect here some of the basic functions used above, as defined in [8]. 
They are the massless 1-loop propagator at zero temperature 

^ 8l ,8 3 ,a)-\p ) J {27T)d[q2]Sl[{q _ p)2]S2 (4n)^T( Sl )T( S2 )T(d-s 12 ) ' 1 ' 
the 1-loop bosonic tadpoles 

^ |Q |» 2T((2s-a-d) T(s - * ) _a 1 . (A 9] 

s % [Q 2 Y (2*T)*-*-* (47r)^r( s ) ' s - % s ' 1 ' ' 

a specific 2-loop tadpole 

Ml, sz,s 3 ; S4 , d) = [Q2]S1[ / 2 ^ 4 _ Q)2]S3 (A.3) 

_ 2T 2 c(2g 123 -2d- S4 ) r(si 3 - |)£(£i2 - |)r(f - si)r(si23 - d) 
^TrT) 2 ^^-^ (47r) d r(s 2 )r(s 3 )r(d/2)r(s 112 3 - d) 

where the shorthand s a foc... = s c + + s c + . . . . 

B IBP relations for 2-loop sum-integrals 

It is known that, for integer indices S1...5, all 2-loop sum-integrals of the form 

T! A\-4 (-Po) S4 (QoY 5 (nu 

L(S1 - 5 ' d) = % Q WFWYWTWF { } 

are trivial in the sense that by systematic use of IBP relations [10] they reduce to products 
of 1-loop sum-integrals (which are in turn known analytically, cf. Eq. (A.2)). We read the 
specific bosonic cases that are needed for the present calculation from our algorithmically 
generated tables, to get 

= - (rf -2Ki-4)(i ( -4-n) ^ W ^ W ■ (R2) 
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L(311,00,d) = ~ {d _ 2){d _ 7) l!(d)I° 2 (d) , (B.3) 
L(221, 00, d) = , (B.4) 

L(211,00,d) = - (d _ 5) 1 (d _ 2) I 2 (rfK2 ^) , (B.6) 

L(211,02,d) = ^|i§(d)J?(d) , (B.7) 
a — 5 

KS11,«,0 = - gzgg^+^ ^CdO . (B.8) 
C IBP relations for 3-loop zero modes 

For the zero-mode sum-integrals 5 denned in Eq. (4.1), for which the generic IBP relation 

= <9 p po 5(si... 5 ;s 6 ,s 7 ) 
= {(d- 2 81 - 8 4 - s 5 ) + s 4 4 + (2^ - 1") + s 5 5 + (3^ - I")} o 5(si... 5 ; a 6 , »?) (CI) 

holds and the additional symmetries <5p II a ft jO dd = and 5p H a b c = (— l) c 5p Ih, ac can be used, 
we give a number of linear relations here that prove useful in Sec. 4 of the main text. Acting 
with the IBP Eq. (C.l) on the sum d ~ 2 2 w+1 5(nllll; ee) + S(n - 1, 2111; ee) (with e an even 
integer) and using symmetries gives the generic relation 

= - (d - 2n + l)(d - 2n - 2)5(nllll; ee) - 5(n - 2, 2121; ee) - S(n - 2, 2211; ee)+ 
+ (d - 2n + l)J|A(nll; e, d) + i|A(n - 1, 21; e, d) , 

which when applied at {n, e} = {3,2} gives 

= - (d - 5)(d - 8)5'(31111; 22) - 5(12121; 22) - 5(12211; 22)+ 
+ (d - 5)lf A(311; 2, d) + if A(221; 2, d) . 

In complete analogy, acting with the IBP Eq. (C.l) on 5(41111), on the sum 5(32111) + 
^ 5(22211), as well as on the sum 5(32222) - (d - 10)5(33212) + ^ ~ 1 9 d + 86 5(33311) and 
using symmetries gives 

= (d - 10)5(41111) - 25(32111) + 2iJM(411; 0, d) , 

= (d - 7)5(32111) - 5(22121) - ^ 5(12221) + i°A(321; 0, d) , 

= - (d - 6)(d - 9)(d - 12)5(33311) + (d - 14)5(23222) - (d 2 - 21d + 106)5(23321) . 
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D Numerical evaluation of finite integrals 



In this Appendix, we treat the finite parts of non-zero (D.l) and zero-modes (D.2) that are 
needed in the main text, cf. Eqs. (3.19) and (4.16). 



D.l Contribution to from finite parts V f of non-zero modes 

Here, we discuss the evaluation of the finite integrals V i as defined in Eq. (2.16), and as 
needed for Eq. (3.19) of the main text. They can be treated at e = in coordinate space. 
Note, however, that we need some of them in shifted dimensions (in fact, for d £ {3,5,7}). 
We will first provide the (inverse, spatial) Fourier transforms of propagators 1 /P 2 as well as 
two-point functions n(P), and then reduce the coordinate-space representation of the V { to 
simple integrals as far as possible, once more in a generic way with unspecified indices. We 
will profit from the fact that we are interested in odd d, in which case the Bessel functions 
coming from the Fourier transforms reduce to Bessel polynomials. The specific cases that we 
need for M3-2 are then approximated numerically. 

The d-dimensional angular averages (Re(d) > 1) result in Bessel functions of first kind 



d-l 



2vr— 



d d rf{r)e ipr = I drr^ 1 f(r) J / du (1 - u 2 ) V e ipru 

Jo r(— ) J-i 

00 

dr r d ~ l f{r) 2{2tt) — (pr) 2 ~ d jdApr) (D.l) 

2 

with j n (x) = ^1 x n J n {x) , (D.2) 
J{I 3 5} = {sin(x), sin(x) — x cos(x), 3 sin(x) — 3x cos(x) — x 2 sin(x)} , 

such that propagators transform into modified Bessel functions of 2nd kind (0<Re(d) <4s+l) 

1 r f A d \c p -* kr r n 2s-d poo kjd ,(k) 

(p 2 + m 2 Y J J (2ir) d (k 2 + m 2 ) s J (2vr)^ Jo {k 2 + m 2 r 2 ) s 

d-2s 

2 l - s 1 /m 2 N 

AW) IT 2 j 

d-2 



(27T) 



with K n {x) = J^-e x K n (x), (D.4) 
f 1 3 3 1 

K±f I 2 6y(x) = < 1, 1 + -, 1 + - + } , (D.5) 



Vr. 



±/a 3 ] 5 ~i — 1 — 2 

zn\ 2 , 2 , 2 f x x x z 
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and we get the Fourier transforms 1 

d d q (Qo) S3 

Qo 



(2ir) d [Q 2 M(P + Q) 2 ]^ 

(1 (r, rp\d+l + S 3 -2si2 



d d re i P r e -\P \r + l^J f J J 2 T Ax (D6) 

j c e -a;(b'l+b'+»»|-|n|) 

with fd,abc( x ' n ) = V— d^i r-53T«: d(|i|x)K 6 _d(|i + n|x) . (D.7) 

i lil a "— + 2 

D.l.l First part V} of Eq. (2.16) 

A generic form for the first part of V i from Eq. (2.16) (plus its {l-H-2} part) is now 

vf(d, Sl ... 8 ) = (n S247 - nf 247 ) (n S358 - nf 358 ) (d.8) 



p (p 2 ) s > 



2^3«-ir( Sl )r( S2 )r( S3 )r( S 4)r( S5 )? J^ x Jo Ve 



-\n\(x+y). 



\4-d-2s! 



x ".,/ WM*. \n\y) j^f T (z,n) / d n -f g (^n) , (D.9) 

2 s r(s) f 00 z 3_d 

hd,s( x >y) = dz TT, — o^^^.^ xz ^d-l(y z ) , (D.10) 

j c e - a: (|j|+b'+™l-l™l) 
a ~~\j+n\ b ~— 

n c y5 r(a)r(6)(47r) ± ? i g c j(a, b, d + 1) K a+b _ d _j_ 1/2 (|n|a;) 
~ M c ^ T(a + 6 - - j) (x/2)i |ri|^+ a + 6 - c - d ' ^ ' J 



In an expansion around odd d and for integer indices Sj, the Bessel functions k reduce to 

■a (a* 
172 



(reverse) Bessel polynomials k s (x) = -^rpw- > such that (a = a + b — d — j — 1/2) 



|o_||_i| 6 _d|_i 



*n-B, , V- (2|a-f|-l-fc)» (2\b - || - 1 - £)! s C!SHa -| ho+ g,, H6 -| h 6+gM 

Jd,abc l X ' n ) ~ 2_> 



*4 d« - II - I - (I 6 - II - 5 - ^ ! (2r E )l a -il+I^H- 1 - fe ^ 

n c r(a)r(6) l^ 1 (47r)^ 5c j(a,6,d+l) 't? (2|a| - 1 - k)\ (|n|x) fc - a -H +c -^' 



|„|c x ,i+c-a-6 ^ T(a + f) ^ (|a| - I -k)\k\ ^\a\-\-k-j 

(D.12) 

where s cab (x, n) = ]T e -*(lil+tf+«|-M) j c \j + n| b = l cafe (coth(x), n) , (D.13) 
j 

Scab(y, -n) = (-l) c s cab (y, n) , s 0ab (y, %) = s ob a(y, %) , s cab (y, 0) = s c>a+bi o(y, 0) , 
s Ccvcn ab(y,n) = s , a +c,b{y,n) , s Coddab (y,n) = s 1>a+c - 1>b (y,n) , ( D - 14 ) 

'From Eqs. (3.9), (3.10), (3.11) we see that /£ oie (x,-n) = M) c f$ abc (x, n) for all X G {II, fl, <7, £>}. 
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or explicitly e.g. 



sooo(c, n) = c + \n\ , 3 io(c, n) = \ [c 2 - I + n 2 + \n\c] , 
s 2o(c,n) = ±[3c(c 2 -l + n 2 ) + |n|(3c 2 + 2n 2 -2)] , 
%i(c,n) = ±[3c(c 2 -l) + |n|(3c 2 + n 2 -4)] , 
sioo(c,n) = -f [c+ |n|] , sn (c,n) = -f [3c 2 + 2n 2 - 2 + \n\3c] 



(D.15) 
(D.16) 
(D.17) 
(D.18) 



Note that g c j(a,b,d + 1)/T( 



a + |) in Eq. (D.12) needs a proper e — )■ limit. Meanwhile, 
the j reduce to trigonometric functions, cf. Eq. (D.2), for which the integral in h evaluates to 
exponentials times polynomials, 



-1) 2 

2e x +v 



0(x-y) [p+ s (x,y) + e 2y p ds (x,y)) + (x^y)+0(e 



(D.19) 



where the polynomials p^ contain even and odd parts of reverse Bessel polynomials 



+ 



k=l 



S - 1 

k - 1 



d-5 (x) 



y a ~ z t 



0,s-k- 



.(y)±y 



2s-2k 



0,k-s+- 



(D.20) 



tb,n( x ) = ( 2n ~ l)*M-i( x ) + x 2 t b , n _ 2 {x) , t fei0 (x) = 1 , = 1 + bx , (D.21) 

h,n<o(x) = X 2n+1 [6 *+,_!(!) + C n _i(x)] , *6,n>o0) = + &*n 0*0 > 



fc=0 



(n-fc)!ife! 2"+! 



■1)' 



(D.22) 



or explicitly e.g. 

pjsfo y) = (y ± !)( 3 + 23; ) ± 0* + x2 + y 2 ) . ( D - 23 ) 

p± 3 (x, y) = (y± 1)(15 + 15a; + 6x 2 + x 3 + y 2 + xy 2 )± 

±y 2 (5 + 5x + 22; 2 ) , (D.24) 
pj 4 (a;, y) = (y ± 1)(105 + 105x + 45x 2 + 10x 3 + x 4 + 10y 2 + lQxy 2 + 3x 2 y 2 )± 

± y 2 (35 + 35x + 15x 2 + 3x 3 + y 2 + xy 2 ) , (D.25) 
pf 3 (x, y) = (y± 1)(315 + 315x + 135x 2 + 30x 3 + 3x 4 + 30y 2 + 30xy 2 + 12x 2 y 2 + 2x 3 y 2 )± 

± y 2 (105 + 105x + 45x 2 + lOx 3 + x 4 + 3y 2 + 3xy 2 + x 2 y 2 ) , (D.26) 



such that (using f(—n) 



-i) c /(«)) 



(47r) 2 



n=l 



1 S / dx / dy 



(47T) 



"^3d+l+ S678 -2s 1234 5, d 



(2ttT) :: 



'vr^-^-l)- 



x [1 + (-1) 5678 ] 



36(xy) 2d ^ 2s 35+^8 



o (4vr 2 ) d -i 2-i2345r( Sl )r( S2 )r( S3 )r( S 4)r( S 5)r(f; 
Pd si i nx ' n y) + e ~ 2ny Pd si ( nx ' x 



j^S6+4— d— 2si g— 2na; 



.r 



2s 3 5-2s 2 4+S7-S g 



+ (x yU + 0(e) 



(D.27) 



-17- 



where fj~f 7 (x,n) = 6x d+s ™ /£f 7 M . (D.28) 

For our sum-integral ^3,-2, we need to evaluate five specific cases which follow from 
Eq. (D.27) and are all of similar structure, such as e.g. (csch 2 (:r) = coth 2 (x) — 1) 

rp2 00 roo i-x g— 2nx . 

4V/(3, 31111, 022) = 7^4 V dx dy ( 3 + 3nx + n 2 x 2 - ny - 2n(l + nx)y+ 

( 47r ) ~\_Jo Jo 18n 5 x 4 y 4 V 

+ n 2 y 2 - e~ 2ny (3 + 3nx + n 2 x 2 + ny + 2n(l + + n 2 y 2 )) x 

x ^3 + nx(3 + 3nx — x 2 ) — 3x 3 (n 2 coth x + (n + cothx)csch 2 x)^ x 

x ^3 + ny(3 + 3ny — y 2 ) — 3y 3 (n 2 cothy + (n + cothy)csch 2 y)^ + 

+ 0{e). (D.29) 

All sums are of the form s a (z > 0) = Y^=i e ~ 2nz I 71 "" = Li a (e _2 ^) resulting in s a> \{z) = 
Li a (e" 22 ), s\(z) = — ln(l — e~ 2z ) and s a <o are polynomial of cothz. For the numerical eval- 
uation of the various Vj we use Mathematica [23] (%2 n — > thousands of LL; then we utilize 
NIntegrate[. . . ,{ x > , 1000} , {y , ,x} , MaxRecursion->100 , WorkingPrecision->20] ). As 
a result, weighting each piece by the prefactor with it contributes to we obtain 



4V/(3, 31111, 022) = 


[+0.01854774(1) + 0(e)] , 


(D.30) 


12V/(5, 31122, 011) = 


[+0.02392697(1) + 0(e)] , 


(D.31) 


9V/(5, 41111, 000) = 


[+0.00006691(1) + 0(e)] , 


(D.32) 


6V/(5, 32111, 000) = 


[+0.00100101(1) + 0(e)] , 


(D.33) 


60V/(7, 33311, 000) = 


^^[+0.01888983(1)+ 0(e)] . 


(D.34) 



D.1.2 Second part V l n of Eq. (2.16) 

The second part of V i as in Eq. (2.16) is (this is only half of it; the other half is 247 358) 
(27rT) M + 3+S678 - 2si2345 T(i) 1 + (-1) 



Vh(d, -1...8) = JC'^T " <*r) «ss ~ n ?3 58 ) (D.35) 



y n d+3+ S68 -2 Sl35 / dxe - 



2nx v 
x 



7T 2 2 2s i234 5 +2d r(§) r(s 2 )r(s 4 ) ^ 

X ^ M " d /Zlfr(^") A-iaBsC^. W' e ) . (D.36) 

N/2] 



f d , si558 (x,y,e)= 2J - " d+1 (l- (^f) ) ± e(4 ? r)+ g S8ij (s 3 , s Bt d+1) 
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In 



xye • 



-d,S1358 ' 



(x) + e 2x Ei(-2x)£ d , Sl358 (-x) + £ d>sl358 (x) + 0(e) , (D.37) 



/A // \ 

(j) = E 9«A*3,*B,d+l) 7 )(d,*136-^ 

w rf , Sl358 > . 'VV 



Q+0(e) 



(<i, s, x) 



X^ 2 ~ c _ d /7T 

e x 2 



r(«) 



„d-2 



2 K l- (x)= r(i) 



poly(x) 

d+i (x) , 



reverse Bessel poly 



n<0 (x) = x 1+2n #_ n _i(x) , e n > (x) = 



£ 2 (d,s,x) = 



"n<0 


(X) 


e x 


POO 




dz 




'o 


In ( 


xe 7 



k=0 



(2n-k)\ (2x) 
(n- ~2» 



(D.38) 



(D.39) 



(D.40) 



In- 



1 + ^ 



(D.41) 



From I and the cases that are useful for us read 



^{3,3112}, {5,3121}, {5,4110}, {5,3210}, {5,3110}, {7,3310} ( x ) = {2x 2 , -2x 3 , -X 3 , 3x 3 , -2x 2 , X 4 }/12 , 
^{3,3112}, {5,3121}, {5,4110}, {5,3210}, {5,3110}, {7,3310} ( x ) = ~~ {2x 2 — 2x, — 2x 3 , —X 3 , 3x 3 , — 4/3x 2 , X 4 }/8 . 

The first part of f u ~ c = f u ~ B - f c ~ B is defined in Eqs. (D.ll) and (D.12) above, while the 
latter part is given by 



i + (-i) c 



d-l d+1 

In 2 x 2 



C(2a- C -rf) ' ' . 2 a r(a-f) 

2 2 7T2 

[c/2] 



(|n|x) 



|n| 6 x a 



-+ 



n=0 



+ [-sign(n)] c E ( 2 C n ) « 2b ~ 2n ~ d )\ n 
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n a _d{\n\x) 



c-2n \ n \ 2 x 2 nb V ( h _ d\ _ 

d-l 1 1 2 J I |a T 6 

2 2 7T2 I ''I x 



We can therefore write the finite parts V\ 1 + 0(e) = Vjj(d, si...g) + Vn(d, S13254687) + C( e ) a s 



4^(3,31111,022) 
12^(5,31122,011) 
9^(5,41111,000) 
6^(5,32111,000) 
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(4vr 
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mvh(7, 33311, 000) = Yl J Q dx y A n 3lo ( 



x, n) 



Ei(-2nx) -e~ 2nx \n \ 
Ei(-2nx) + e~ 2nx In 
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Note that line 3 and 5 cancel exactly (assuming they are finite). The reason is that n§ l0 = 
v^ 2 ' 1 ^ = 3^1 goo(i,M+i) = 2^nf in such that the relevant contributions to M 3 -2 

[p2]3--3- f [P2] 2 — 5- ^ 

combine with a pre-factor of 0(e): 



M™f 2 3 2dV + |3V/ 7 (41111, 000, o 3 ) + V/^31211, 000, o 3 )} 



= 2dV + 13^(41111,000,03) + (2 - d) 1^(41111, 000, a 3 )} 
= 2dV + j (5 - d)F/ 7 (41111,000,a 3 )} 

= 2(f(3 - d)V+Vfj (41 11 1,000,03) . (D.43) 

Thus, also the contributions from V d should not contain 03 up to their constant terms, which 
is indeed the case, serving as a small check of our expressions. Furthermore, it appears 
convenient to choose oi = 02 = 05 = 16tt 2 e 3 / 2 " 7 and 04 = 16n 2 e 1-7 , although the 7 could 
also remain in the log, see Eq. (D.52). 

The various functions f u ~ c above are (omitting the arguments (x, n) on the lhs) 

?n-C_ 1 (i 1K . k / 2 on K 2 2/o , 2x , -,^3/2 



/£~2 = __ ^ 4 _ 15 + 5 nx ( x 2 - 3) - 5^x^(3 + x l ) + Vox 6 (n z coth x + (n + cothx)csclTx 



/ipm = ~^~~2 {p + + n(3 + x 2 )) — 3xcothx(l + nx + xcoth 
ff-g = (225 + 135nx + 15nx 3 + x 4 + nx 5 - 

— 45x(2(l + nx) cothx + x(2 + nx + x coth x)csch 2 x 

= -— Ij^ ( 270 + 120x2 + x4 + 30nx ( 3 + - 90x coth x (2 + nx + x coth x)) , 
= - 1 3 (9450 + x(3780x + x 5 + 315n 2 x(3 + x 2 ) + 315?i(12 + 5x 2 ))- 

— 945xcothx(6 + 3nx + (—1 + n 2 )x 2 + xcothx(3 + nx + xcothx)) 

A quick numerical evaluation (summing up to n max = 10000, and using Mathemat- 
ical NIntegrate[. . . , ,MaxRecursion->100,WorkingPrecision->60,AccuracyGoal->30] , 

while estimating the error by fitting values n G [10000, 60000] and extrapolating to infinity) 
of the (weighted) Vjj parts then produces (choosing all Oj = 167r 2 e 3//2_7 ) 

T 2 

4^(3,31111,022) = —-j [+0.00775440(1) + 0(e)] , (D.44) 
(47r) 4 

T 2 

121/^(5,31122,011) = j—y [-0.00354681(1) + 0(e)] , (D.45) 
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rp2 

9V}j(5, 41111, 000) + 6V/ 7 (5, 32111, 000) = 7—= [+0.00295006(1) + 0(e)] , (D.46) 

(47TJ 

T 2 

60V/ 7 (7, 33311, 000) = Q [-0.00503877(1) + 0(e)] . (D.47) 

It might be possible to evaluate (some of) the Vjj analytically. We have not put further effort 
into this, as the numerical values given above are fully sufficient for our purposes. Before 
leaving, let us record some integrals and some sums that might be useful in that respect: 

f'°° f'°° F(n + 1) 

dze~ z z n = T(n + l) , /dzEi(-z)z n = tSlLL±l , R e (n) > -1 , (D.48) 
Jo Jo n + 1 

/•oo 

/ dz e~ z z n (7 + In z) = T(n + 1) HarmonicNumber(n) , Re(re) > -1 , (D.49) 
Jo 

r°° 1 7r 2 

dz [e- z in + Inz) - Ei(-z)] - = C(2) = — , (D.50) 
JO Z 

E^ = -C», E^=Li a (e--), t e -^^ = -O a U a (e^). (D.51) 

n=l n=l n=l 

One can get some analytic parts from looking at the piece ~ fsno (x,n) = <|j (1 + \n\x) of 

T 2 f f? fs,noM 



^(5,41111,000,03) dx 



18n 



7—^ — Y \ — — 7 — /dz z(2 + z) [e" z ln(z e 7 ) - Ei(-z)l - 
(4vr) 7 18 ^ \ 108n 3 20 7 V ' L V ; V 7J 

ra-1 s y / 

=1 

Inn 1 , s z ) T 2 1 /C(3) C'(3)\ 



90n 3 4 



where in the second line we have chosen = 167r 2 e 3//2 , and which is confirmed by the 
corresponding analytic expression for the (U c - U B ) piece of V d , which can be extracted by 
considering only the I ■ I parts of Eq. (3.15). 

D.1.3 Summing up 

Summing up Eqs. (D.30)-(D.34) and Eqs. (D.44)-(D.47), we obtain the numerical coefficient 
of Eq. (3.19) as 

m » +0.0645513(1) . (D.53) 

D.2 Contribution to from finite parts V z ' of zero-modes 

Here, we discuss the evaluation of the finite terms given in Eq. (4.16) of the main text. 
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Working again in coordinate space, we need the Fourier transform for Tl A of Eq. (4.8). 
Being the sum of two terms, its first part can be read off from the last line of Eq. (D.12), 
while in the notation of Eq. (D.6) with Pq = its second part (let us label it E here) reads 



Y(d/2 - a)T(d/2 - b) /2vrTr\ 



a+b+l-d 



1-2(1/2) 



— 



(D.54) 



whereas below we also need the special case /•f 2 22( x ' 0) = 1- 

What we need to evaluate is the first line of Eq. (4.7) for which, using Eq. (D.2) we obtain 



Sp 

[p2]«i 
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r( Sl )r( S2 )r( S3 )r( S 4)r( S5 )r(d/2)^i+3 d / 2 J Q dx J*v h ^ (*> v) x 



„S24" 



y 



with h d>3 (x, y) = r( ^ 2 r (d/2) d 2 ^ 2Fl{d/2 " s ' 1 " s ' d / 2 ' y2 / x ^ 
and special cases ^{31,51,52,72} (x, y) = {y, ^, ^ 5x ^ y2 , I^ 3 ^} . 
The function h in fact originates from 



(D.56) 
(d>2s>0) , 

(D.57) 



2 s r(s) 



7T 



dzz 3 d 2s j d /2-i(x z)j d/2 -i(y z) = 6(x - y)h djS (x,y) + 9(y - x)h d)S (y,x) 



which corresponds to the overall J of Eq. (D.55) after performing the angular integrals. 

To get an explicit expression of f^~ A , we simply refer to Eq. (D.12), taken at n = and 
where we only have to insert the theta function of Eq. (4.8) into the second line, as well as 
Eq. (D.54) above, again multiplied by the theta function. For the specific cases that we are 
interested in, we get the compact expressions 
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f 5ioro(3j 12211, 22 



-f SUA 22121, 00) 



fi(x,y) = x 5 (cothx — 1) (y 3 coth(y)csch 2 y — l) , 

j>2 



poo px 2 

dx dy — 7T- (l — x 2 csch 2 x) (l — y 2 csch 2 y) (D.59) 
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h(x,y) = x 3 (cothx - 1)( - 3 + ycsch 2 (y) (y + sinh(2y))) , 
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/sOej y) = x 3 (coth x — 1) (2+x+x coth x) (6y coth y — 10+y 2 (3+y coth y)csch 2 y) 
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fe(x,y) = x 3 (cothx — l)(2+x+x coth x) (2y cothy — 5+y 2 (2+y cothy)csch 2 y) . 

The y-integration can be explicitly performed in most (all but the last two) cases, giving 
a result containing zetas and (poly) logarithms, whereas the remaining x-integration can be 
approximated numerically, to give 



|5 f (3, 12121,22) p 
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" (4tt) 4 


[ - 0.01005114745(1) + 0(e)] , 
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[ + 0.15213227462(1) + 0(e)] , 
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" (4tt)^ 


[ - 0.00891125885(1) + 0(e)] , 
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-^5 f (5, 12221,00) p 
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[ - 0.10673211253(1) + 0(e)] , 
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96 ^(7, 23321, 00) p 


S (^FT [ + 0.20021689747(1) + 0(e)] , 
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84^(7,23222,00) p 


* (^F [ + 0.02318282360(1) + 0(e)] . 
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D.2.1 Summing up 

Summing up Eqs. (D.64)-(D.69) we obtain the numerical coefficient of Eq. (4.16) as 

n 2 « +0.24983747686(1) . 
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E Cross-checks 

Using the generic formulae for spectacles-type sum-integrals as derived in the main text, we 
can check two particular cases, Mi^ and V2, against the previously known results from the 
literature. 

First, let us check Aitfi of Sec. 2 in [11]: In the language developed here, 

Mi,o = 7(11111, 000) = V 1 + V d + V z = B + 2C + y d (lllll, 000) + A + S d (lllll, 00) , 

with {A,B,C} as defined in [11] and whose expansion (using p 31 (x, y) = ±1 as well as 
Eqs. (3.14), (4.9) and £(111, 00, d) = 0) exactly reproduces Eqs. (2.15)-(2.17) of Ref . [11]. 
Second, let us check V 2 of [12]: In the language developed here, 

V 2 = 7(21111, 002) = V 1 + V d + V z = V 2 + V 4 + V 6 + y d (21111, 002) + 5(21111, 02) , 

where IBP on 5(21111,02) gives the reduction 

Vi = 5(21111,02) = - — - (5(12111,02) + 5(12111,20) -i§A(211,2) - if 4(211,0)} 
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for which 



5(12111,02) = Via + 5 d (12111, 02) , 
5(12111,20) = Vi 6 + 5 d (12111,20) . 



(E.l) 
(E.2) 



Collecting and expanding (using p 32 (x,y) = y ± (x + 1); Eqs. (3.14), (4.9); as well as 
L(lll, 00, cf) = 0; L(211,20,d) = 0; and L(211, 02, cf) from Eq. (B.7)), one reproduces 
Eqs. (4.1)-(4.2) of Ref. [12]. 
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